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Correlations and Fluctuations, A Summary of Quark Matter 2002 
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Results for correlations and fluctuations presented at Quark Matter 2002 are summa- 
rized. These results include Hanbury-Brown Twiss interferometry of a wide variety of 
species, large scale fluctuations and correlations in pt and multiplicity, and charge fluctu- 
ations and charge balance functions. 

1. Introduction 

A phenomenal number of new results have been presented at this meeting regarding 
fluctuations and correlations. With the analysis of SPS results for tt — vr Hanbury-Brown 
Twiss (HBT) correlations as a function of energy, and with the maturing of the RHIC 
results, a consistent picture has developed of the evolving dynamics from AGS energies 
to RHIC. Additionally, several new probes have been analyzed and have been presented 
at this meeting, including 77 correlations. Fluctuations in pt, multiplicity and isospin, 
derived from both the SPS and RHIC, hint at novel behaviors. Finally, charge correla- 
tions measured at RHIC suggest a delayed production of charge, and therefore a delayed 
hadronization. After being taken by surprise, the theory community has made signifi- 
cant progress toward understanding these measurements and reducing the list of possible 
interpretations. 

Results on small relative momentum correlations are reviewed in the next section, 
with an emphasis on describing the HBT puzzle. The subsequent section reviews pt 
and multiplicity fluctuations, while isospin fluctuations, charge fluctuations and charge 
balance functions are reviewed in Sec. lU 



2. Two-Particle Correlations at Small Relative Momentum 

A principal motivation for HBT analyses was the hope to view a long-lived source 
resulting from the low pressure associated with a flrst-order phase transition to the the 
quark-gluon plasma (QGP) [[I1I21E]- The CERES collaboration presented results, shown 
in Fig. ^ of both their results at the SPS along with AGS and RHIC results [IH El- A 
long-lived source would have emitted pions continuously over a time of the order of 10-20 
fm/c, which would result in the outward size i?out being larger than -Rside by an amount 
of the order of 10 fm. Figure [H shows that such a long-hved source was not produced at 
any energy. In fact, the outward and sidewards dimensions show little dependence with 
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Figure 1. The three dimensions of the source are illustrated as a function of transverse 
momentum for several beam energies. Only -Riong changes significantly for higher beam 
energies. 



beam energy. Consistent measurements of the three source dimensions at RHIC were 
reported in this meeting by STAR, PHENIX and PHOBOS, with the ratio of -RoutZ-Rside 
varying from 0.85 for PHENIX [|^ to 1.15 for PHOBOS [Ej. Because of the inherent 
ambiguity between the temporal and spatial dimensions in -Rout, one can only place an 
upper limit on the duration of the emission of Ar < 10 fm/c. Given the large transverse 
size, it appears that the breakup hypersurface is extremely space-like, i.e. emissions from 
the edge and the center occur within a few fm/c of one another. 

The values reported for -Riong are as remarkable as the i?out/-Rside measurements. Whereas 
i?out/-Rside provides insight into the duration of emission, -Riong gives a measure of the emis- 
sion time. HBT measures the spatial extent from which pions of a given rapidity might 
be emitted. Assuming that the emitting sources are well spread out in rapidity, the size 
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Figure 2. Theoretical predictions of Kolb et al. compared to STAR results. As indicative 
of several other treatments, the inclusion of the final hadronic rescattering pushes the 
theoretical predictions further from the results. 



is determined by the temperature and the velocity gradient, -Riong ~ v^h/ {dv / dz), where 
Vth is the thermal velocity and dv/dz is the velocity gradient. In a Bjorken expansion, the 
matter does not accelerate along the z axis which means that the velocity gradient is de- 
termined by r, dv/dz = 1/r. Although the values for -Riong reported by STAR, PHENIX 
and PHOBOS are very large compared to hadronic length scales, they still fall short of 
expectations, and suggest emission times of 10 fm/c, rather than the expected 15-20 fm/c. 
Unless the breakup temperatures are anomalously cold, or unless the matter accelerated 
in the longitudinal direction, explaining this result requires an extremely rapid transverse 
acceleration. Blast wave analyses [ |8j of tt — vr HBT reveal sources which have expanded 
to an outer radius of 13 fm within this 10 fm/c, and with an outer velocity of ~ 0.7c as 
determined by comparing pion and proton spectra. Given that the original radius is near 
6 fm, this implies an extremely rapid acceleration. 

The inability of numerical models of the dynamics to account for both the small values 
of the mean emission time, (r), and the duration of the emission. At, is referred to as the 
HBT puzzle. Several attempts have been made at reproducing the RHIC results using a 
variety of hydro dynamic, Boltzmann and hybrid treatments. All consistently overestimate 
i?iong and -RoutZ-Rside [El UHl UH 112] • Figure |21 displays results from Kolb which come 
somewhat closer to the result than the other reference work cited above. It should be 
emphasized that the HBT puzzle is not particularly new to RHIC. As HBT statistics 
have increased form the SPS, it is clear that the emission duration was moderately over- 
predicted by RQMD at SPS energies, with measurements of -Rout/-Rside varying from near 
unity for CERES [ij to ~ 1.15 for NA49 [HH]. 

Two-kaon correlations from PHENIX and from NA49 were also reported at this meet- 
ing [ El . It appears that the source size for kaons is somewhat larger than predicted 
by extrapolating from the pion source sizes assuming blast wave kinematics. However, it 
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should be pointed out that the NA49 measurement is somewhat larger than that previ- 
ously reported by NA44 [E], and that the PHENIX measurement is preliminary. 

Aside from the kaon measurements, the overall behavior of HBT measurements appears 
to follow the systematics of a blast wave, i.e. a sudden dissolution at a low temperature, 
~ llOMeV, after a rapid expansion to a high velocity, ~ 0.7c. In addition to the -RoutZ-Rside 
ratio and the small value of -Riong, the blast wave picture of sudden explosion and emission 
explains the pt dependence of the source sizes [ [UJ [T21 IH CHI, and may also explain the 
results of STAR where two-pion HBT was performed relative to the reaction plane [ 
Unj. There, the source appears to be out-of-plane extended, i.e. the source appears to 
retain the shape of the original almond resulting from overlapping spheres at non-zero 
impact parameter. If the emission were long-lived, the differential expansion would have 
overcompensated for the initial deformation and resulted in an in-plane-extended shape. 

One consequence of the smaller-than-expected values of -Riong and i?out is that the phase 
space density must have been rather high at breakup. In fact, by combining the spectra 
and HBT measurements, one can extract the phase space density in a model-independent 
way [IF]. This analysis was performed for STAR results [El, and the resulting phase 
space density suggests that pionic phase space is overpopulated by nearly a factor of two 
at low momentum as compared to a thermally equilibrated gas where the pion density 
is determined solely by the temperature. The phase space density increases markedly 
with centrality and is significantly higher than was extracted for central collisions at SPS 
energies [El. It should be emphasized that the reported phase space values are averaged 
over coordinate space, 

f( N ^ Jd^xfjx, p)/(x, p) 
^^P^- /rf%/(x,p) ' 

For a Gaussian source the peak phase space density is larger by 2^/2 at x = than the 
average derived from the sum rule. This suggests that the break up density is a significant 
fraction of what is required for novel Bose coherence effects. 

A high phase-space density necessitates a lower entropy. In a nearly model-independent 
fashion, the entropy carried by the pions was extracted from STAR results [ E] • The 
entropy per pion appears to be in the neighborhood of three units which is significantly 
lower than that expected for a Bose gas as can be seen in Fig. 01 As hydrodynamic 
treatments over-predict the source sizes, while fitting the spectra [ E] , they thus over- 
predict the entropy. Since hydrodynamic treatments produce minimal entropy in their 
evolution, one must conclude that the initial entropy of the system in the hydrodynamic 
treatments was overestimated. This suggests that the initial state was not that of an 
equilibrated quark-gluon plasma, but might have instead involved energy in coherent 
forms, e.g. classical gluon fields, or may have excited fewer degrees of freedom, e.g. 
under-populating quark degrees of freedom. 

The correlation functions presented at this conference for non-identical particles were 
shockingly superior to what had been previously measured. Following the ideas of Led- 
nicky [ 20j, correlations were compared for p^^out — Pfe.out > to those with p^^out — Pfe.out < 0. 
Whereas the two correlations must be identical when the particles are indistinguishable, 
the difference offers insight into the relative positions of the emitted particles when the 
two particles are different. STAR presented results for tt — i^T and n — p pairs [|H1 which 
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Figure 3. The entropy per pion as extracted from HBT and spectra. For the most 
central colhsions, the entropy per pion appears low which puts doubt into assumptions 
about early equilibration. 



showed that the heavier particles were emitted ahead of the lighter particles consistent 
with a blast wave interpretation. Results from NA49 for tt — p pairs demonstrated the 
same behavior [IT^. 

In this meeting correlations were reported for several new pairs which had not previously 
been measured in heavy ion collisions. NA49 displayed results for both AA and Ap 
correlations [E]. The AA results showed a flat correlation, which means that they are not 
useful for source-size determination, but can put limits on the AA scattering length. This 
is important as it sheds light on whether there is a low-energy resonance associated with 
the Ho di-baryon. The A — p correlation function exhibited a peak near zero momentum 
which was consistent with expectations for a three to four Fermi source, though the size 
was difficult to ascertain without knowing the A parameter which depends on the fraction 
of As from decays of heavier hyperons. With greater statistics, this measurement offers 
the possibility of determining whether As are emitted ahead or behind the protons. 

The Coulomb interaction between the residual source and the pairs used for HBT 
has always clouded the interpretation of experiments. However, Kg — Kg interferometry 
sidesteps the issue since the Kg is neutral. A first measurement of the Kg — Kg correlation 
function was presented by STAR [CHI- Although the statistics were insufficient for a 
detailed analysis, this preliminary result was extremely encouraging. 

The most surprising result of the conference was the presentation of two-photon HBT 
by WA98 [ El 122] • The height of the correlation function shown in Fig. |3] has a height 
which requires that 10% of the photons must be direct photons, i.e. not from ttq decay or 
dalitz decays from long-lived particles such as the rj. In fact, by measuring the inclusive 
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Figure 4. The two-photon correlation function from WA98 is shown with and without 
cuts for cluster splitting in the calorimeter. 



photon spectrum, one can use the intercept of the correlation function, 1 + A, to infer the 
direct photon spectrum, 

^ = (2) 

This result was surprising because the A parameter was expected to be much smaller. If 
the fraction of direct photons to inclusive photons were 1% rather than 10% the height 
of the correlation function would have been smaller by a factor of 100, and the statistical 
requirement for events would have increased by a factor of 10^. 

I conclude this section with a brief summary of theoretical attempts at solving the 
HBT puzzle. Most importantly, it should be emphasized that there is consensus that 
hydrodynamics with a Boltzmann or cascade afterburner will over-predict source sizes [ 
El HD 1121 Cni I21I21- However, it is not yet clear that fitting the data would not be 
possible if an extreme equation of state were used, e.g. one with a speed of sound greater 
than 1/a/3. Several careful studies have been performed to ensure that experimental 
resolution and Coulomb corrections have been properly treated [HJ 123 12^1 • The effects 
of viscosity have just begun to be included, and it seems that they will provide only a 
modest correction in order to be consistent with elliptic flow results which preclude large 
viscosities [123 EH]- The effect of modifying chemical rates and cross sections has also 
been investigated with only modest changes in HBT results [ |2S1 ^] ■ Finally, the efforts 
at comparing theory to experiment have become more meaningful as theoretical models 
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have generated correlation functions which are then fit for Gaussian parameters, rather 
than comparing radii generated from looking at source functions. It appears that the 
HBT puzzle indeed derives from comparing "apples to apples". 

In addition to listing all the investigations presented in this conference, it is worthwhile 
to list the issues not yet addressed which might shed light on the HBT puzzle. First, it 
should be pointed out that the hydrodynamic treatments all assume a Bjorken space-time 
geometry, i.e. the matter coasts along the beam axis. If the matter accelerates significantly 
due to the finite size in the beam direction, the extrapolations for (r) at breakup might 
be too small. This issue could be settled by three-dimensional hydrodynamic studies. 
Secondly, three-body interactions such as the residual Coulomb interaction have not been 
totally ruled out as being unimportant, although schematic calculations suggest they 
are indeed unimportant, and the equivalence of vr+vr"'" and 7r~7r~ results also points to 
their insignificance. The most interesting possible explanation of the -RoutZ-Rside puzzle 
might lay in a reduced emissivity of the surface, which would allow the matter to fall 
apart suddenly. The emissivity might be reduced in the case of a super-cooled phase. 
The problem of how quarks might escape microscopically from a QGP region and form 
mesons has been considered from the point of fissioning flux tubes [|30j. Effects of an 
increased emissivity were reported on at this conference [ 122] ■ As a final possibility, it 
may be that the picture of the dynamics assumed in cascade or Boltzmann algorithms is 
fundamentally flawed due to the neglect of quantum effects. Again, this seems unlikely, 
but it should not be forgotten that the thermal wavelengths at breakup correspond to 
volumes with several particles. In order to explain the HBT puzzle, these issues must 
somehow account for corrections of the order of 30%. 

3. Multiplicity correlations and pt fluctuations 

Fluctuations in multiplicity or transverse momentum have been proposed as signatures 
for phase separation [|M1E2] or critical fluctuations [E^. Of course, such fluctuations are 
also affected by jets and jet quenching. In this conference we have seen new results for 
Pt fluctuations from NA49 [US], PHENIX [OH and STAR The results from NA49, 

shown in Fig. illustrate that the fluctuation is smaller for same-sign particles. This may 
be due to the positive correlation between opposite sign particles from resonance decays. 
HBT effects, jets and the distorted shapes of the nuclei are other potential correlations 
which would give pt fluctuations. 

The value of increases if there is an increased correlation in pt between given particles 
or if the number of particles with which one is correlated increases. By measuring the 
kinematic range of the correlation, one should gain insight into the source of the non-zero 
$pj. The fluctuations measured by PHENIX were even smaller than those reported by 
NA49, while the fluctuations measured by STAR were significantly larger. There seems 
to be a trend that the fluctuations are larger for larger acceptances, which suggests that 
the correlation may span more than a unit of rapidity. 

A wavelet analysis, which can measure the same class of correlations as pt fluctuations, 
was presented by the STAR collaboration in this meeting [EH]- For low multiplicities, 
the correlations were consistent with HIJING with jets, while at high multiplicities, the 
results were smaller and were consistent with HIJING without jets. This emphasizes the 
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Figure 5. Fluctuations in pt from the SPS are shown as a function of centrahty for all 
charges or for same-sign particles. 



point that jet quenching should be manifest in low pt observables. Perhaps, experiment 
could decide not only the existence but the nature of jet quenching, i.e. determining 
whether the jet energy is transformed into more particles through radiation or whether it 
is shared among the collision partners. 

4. Isospin fluctuations and correlations arising from charge conservation 

Large fluctuations in the fraction of neutral pions to total pions has been suggested 
as a signal for coherent emission of pions which might result from a disoriented chiral 
condensate [EHl- Unfortunately, this is a tremendously difficult measurement due to the 
fact that neutral pions decay into two photons, making the exact counting of soft pions 
extremely difficult, especially in a central collision. WA98 presented an analysis showing 
no evidence for large fluctuations [|^. However, one should not jump to the conclusion 
that the observations were entirely consistent with random emissions until a more detailed 
theoretical investigation is performed. In this conference, new ideas were also presented 
for investigating fluctuations of a different flavor, i.e. fluctuations in flavor [ and 
baryon number [I38j. 

Charge fluctuations, charge correlations and balance functions represent another class 
of observables, which are sensitive to the dynamics of separating conserved charges. Since 
charge conservation is local in space-time, each positive charge is accompanied by a nega- 
tive charge. If these balancing charges remain close to one another in coordinate space at 
the time of their emission, they will be correlated in momentum space. This means that 
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Balance Function, AutAu, ^ij^ = 200 GeV, Charged Pairs | All Charged PalTs^ 




Figure 6. Charge balance functions (left panel) measured by STAR show that balancing 
charges are more correlated in central collisions than in peripheral collisions. The widths 
are plotted as a function of centrality in the right panel. This behavior is consistent with 
a delayed hadronization. 



the net charge measured into a large volume of phase space, e.g. the STAR acceptance, 
would tend more toward zero since each charge would be canceled by a balancing partner. 
A delayed hadronization would lead to more highly correlated balancing charges since the 
charges would have less time to separate and would also be produced after the early stage 
of the collision where there is a high velocity gradient pulling particles apart. Since most 
charge is created at hadronization, a tighter correlation of balancing charges could serve 
gnal of QGP formation [EHlllI!. 
Several analyses were presented at this meeting which measured this correlation. Charge 
fluctuations were measured by NA49 [E], and charge fluctuations and balance functions 
were presented by STAR The balance functions from STAR are shown in Fig. IHl The 
balance function statistically measures the distribution of relative rapidities of balancing 
7r+ — TT^ pairs. The narrowing of the balance function as a function of centrality is highly 
suggestive of a delayed hadronization. 
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